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Abstract: We reconsider the realization of baryons in the Sakai-Sugimoto model. In this 
theory, which is the gravity dual of a QCD-like theory, baryons appear as soliton solutions. 
These solitons were approximated as flat-space instantons in previous studies. However, 
with this approximation, it has been shown that one does not reproduce some model- 
independent predictions for the behavior of baryon electromagnetic form factors which are 
connected with long-range pion physics. This made it appear that the long-range pion 
physics of baryons may be hidden in (intractable) a' corrections in the gravity dual. In 
this paper, we study the long-range properties of baryons in the Sakai-Sugimoto model 
without relying on the flat-space instanton approximation. The solution we obtain gives 
the correct result for the model-independent ratio of form factors, implying that the model 
captures the expected infrared properties of baryons without the need to go beyond the 
leading order in the a' expansion. 
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1 Introduction 

Thanks to the development of gauge-gravity duality [1-3], we can use gravity duals as 
powerful computational tools for analyzing strongly-coupled QCD-like theories. Perhaps 
the most prominent QCD-like theory with a gravity dual is the Sakai-Sugimoto model[4, 5]. 
The model has a spectrum of weakly-interacting mesons with the same quantum numbers 
as in QCD, and baryons can be thought of as solitons of the meson fields, just as one 
would expect in large Nc QCD [6, 7]. The price paid for the tractability of the model 
is that in addition to the fields and particles that are contained in large Nc QCD, there 
are also many other particles at the same mass scale as the particles with QCD quantum 
numbers. These extra fields conspire to make the 't Hooft coupling A = Qym^c a tunable 
control parameter of the model, which is necessary for the existence of a classical gravity 
dual. Since the Sakai-Sugimoto model contains a great many of the ingredients of QCD, 
it is hoped that for many observables the model will give at least qualitatively accurate 
predictions. Of course, this begs the questions of which observables in QCD have a good 
match in the Sakai-Sugimoto model, and how to define what counts as a good match. 

At least one setting where such questions should have sharp answers is for low-energy 
QCD observables, where there are predictions from chiral perturbation theory (x-PT). 
For such observables it is known precisely what one should expect from QCD, and sharp 
comparisons with the Sakai-Sugimoto model are possible. The question we focus on in this 
paper is whether the predictions of X"PT for baryons are satisfied by the Sakai-Sugimoto 
model. Since the Sakai-Sugimoto model has spontaneous chiral symmetry breaking, and x- 
PT is simply a systematic way to work out the the consequences of symmetry breaking for 
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low-energy observables, this may seem like a foregone conclusion. Things are not so simple, 
however. As an effective field theory (EFT), X"PT assumes that the low-energy constants 
appearing in the derivative expansion are of natural sizes. The Sakai-Sugimoto contains the 
extra parameter A compared to QCD, and it could be that some low-energy constants are 
suppressed by powers of A. If this happens then there will be low-energy observables that 
behave qualitatively differently between QCD and the Sakai-Sugimoto model. Fortunately, 
no observables that disagree with the predictions of X"PT have been found in the meson 
sector. But in the baryon sector^, things are less encouraging, and observables which appear 
to be afflicted by this issue have been previously identified in Ref. [15]. The observables 
discussed in Ref. [15] are the long-distance limits of the Fourier transforms of the electric 
and magnetic form factors of the proton, which encode the response of the proton to 
soft photon probes^. The leading long-distance behavior of these observables is tightly 
constrained by X"PT, and so must satisfy some model-independent relations in QCD-like 
theories. 

Specifically, Ref. [15] argued that in the chiral limit = 0, these position-space form 
factors must obey the relation^ 

where the Cs are isoscalar/isovector electric/magnetic form factors of the proton in posi- 
tion space. These observables are just the Fourier transforms of the usual momentum-space 
form factors and are defined as 

Gr°(0 = ^ / dn{pt\JLo\p^) (1-2) 

Gi7°(0 = i- I dn^s,,3{p t \x.JLo\p t) (1-3) 

G'E=Hr) = 1^1 dn{p t I Jtt'"='lP t) (1.4) 

Gifir) = i- 1 dn^-e,,s{p t |x. Jn'"='lf t) (1-5) 

where Jj^q, Jjtfi the isoscalar and isovector currents, and the matrix elements are 
with respect to a proton in a spin 'up' state to be specific. In the chiral limit the pions are 
massless, and so give the leading contribution to the long-distance properties of the proton. 
The long-range behavior of these form factors is determined by the hadronic diagrams in 
Fig. 1. Note that these diagrams are not suppressed in the large- A'^c limit despite seeming 
to contain a pion loop, because meson-baryon coupling constants scale with [17-19]. 



^For studies of baryons in the Sakai-Sugimoto model, see e. g. Refs. [8-14]. 

^Unlike other probes of the long-distance behavior of baryons such as charge radii, which are related to 
the moments of form factors, the Fourier transforms of the form factors remain finite in the chiral limit. 

^The chiral limit and the large Nc limit do not commute. In this paper we assume that the large Nc 
limit is taken first. While we work at m,r = throughout, it would also be interesting to explore what 
happens at finite m^, as was done in the Pomarol-Wulzer model in Ref. [16]. 
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(a) Isoscalar 



(b) Isovector 



Figure 1. Some representative diagrams that contribute to the long-distance behavior of form 
factors in large Nc X-PT. Fig. 1(a) depicts a leading contribution to the nucleon isoscalar form 
factor, which is dominated by a three-pion interaction from the anomaly. Fig. 1(b) depicts a leading 
contribution to the isovector nucleon form factor, which is dominated by a two-pion interaction. 
Since the A baryon and the nucleons are degenerate at large the intermediate states can be 
either nucleons or As. The As are shown as double lines above. 

The relation in Eq. (1.1) is known to be satisfied by any chiral soliton model of baryons, 
such as the Skyrme model and its many generalizations. It is also satisfied by the Pomarol- 
Wulzer bottom-up holographic model of QCD[20-22]. 

The surprising conclusion of Ref. [15], however, was that baryons in the Sakai-Sugimoto 
model did not appear to obey the model-independent relation. This conclusion was based 
on the standard treatment of baryons as holographic flat-space instantons in the Sakai- 
Sugimoto model [9, 13]. The flat-space instanton approximation is motivated by the large 
A limit, and appears to imply that the pion cloud around the nucleon is absent at large 
A. Instead, the flat-space instanton approximation suggests that the leading long-distance 
contribution to the form factors of the proton in the Sakai-Sugimoto model is mediated by 
p mesons. As a result, the long-distance behavior of the form factors is very different com- 
pared to what would expect from X"PT. The large A limit is mapped to the a' expansions 
in the dual string theory, and as was discussed in Ref. [15], the physics of the pion- nucleon 
interactions is then apparently hidden in the a' corrections in the string theory, which are 
not a calculable part of the model. These considerations suggest that the large distance and 
the large A limits do not commute, so that if the large A limit is taken first the long-distance 
physics is completely different than in QCD. So if the standard treatment in the literature 
is correct, then for this class of observables there are quite sharp qualitative differences 
between the Sakai-Sugimoto model and QCD. Our aim is to investigate whether this is 
indeed the case. 

In this paper, we reexamine the calculations of the long-distance properties of baryons 
in the Sakai-Sugimoto model. In contrast to previous treatments of the model, we do 
not start by assuming the flat-space instanton approximation, and directly solve for the 
long-distance behavior of the holographic soliton associated with the baryon in the dual 
field theory. We then show that the structure of this long-distance solution is such that 
Eq. (1.1) is in fact satisfied by the baryons in the Sakai-Sugimoto model. Our results are 
obtained directly from the leading terms in the action in the a' expansion. So contrary 
to the pessimistic story above, the long-distance pion-nucleon physics is not hidden in the 
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01 corrections to the Sakai-Sugimoto model. It is present in the Sakai-Sugimoto model at 
leading order in the a' expansion. 

This paper is organized as follows. In Sec. 2 we give an overview of the Sakai-Sugimoto 
model, describing the dual gravitational theory. In this overview we emphasize the impor- 
tance of introducing a UV cutoff in the gravity dual, which is not always appreciated. 
We also describe the self-consistent ansatz for the soliton fields which is appropriate for 
investigating nucleons. We solve the equations of motion of the soliton in a 1/r expansion 
in Sec. 3. These solutions are turned into predictions for the long-distance behavior of the 
form factors using the standard machinery of collective coordinate quantization. The form 
factors are seen to obey Eq. (1.1). In Sec. 4, we make a few brief remarks about some 
possible issues in the previous treatments of baryons in the literature, which may help 
explain why the results we obtain in this paper were not obtained in previous analyses. 
Finally, we outline some possible directions for future work. 

2 The Sakai-Sugimoto model 
2.1 Brane construction 

The Sakai-Sugimoto model describes the strong-coupling physics associated with a system 
of intersecting Nc -D4 branes and N ^ D8, D8 branes in type IIA string theory. It is assumed 
that Nf <C Nc- At weak coupling, an open string picture is appropriate, and gluons can be 
seen to arise as open strings with both ends on the Nc DA branes, while quarks arise as open 
strings with one end on the Nc color branes and another end on the Nf flavor branes. The 
direction is chosen to be a circle, with anti-periodic boundary conditions for fermions, 
which breaks supersymmetry. The 4D field theory lives on the intersection of the world 
volumes of the color D4 branes and the flavor D8 branes. Finally, baryon operators appear 
as DA branes wrapping the xg, • • . , xg directions. The brane configuration with Nb baryons 
is summarized in Table 2.1. 
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(2.1) 



We will be interested only in configurations with Nb = 1 in this paper. Ref. [4, 5] showed 
that despite the breaking of supersymmetry by the boundary conditions on the S^, this 
brane configuration is stable, at least to leading order in a Nj/Nc expansion. 

Let us briefly review the well-known connection of this brane construction to QCD[4, 
5, 23]. Because of the supersymmetry-breaking boundary conditions on S^, all fields except 
the gluons and quarks pick up masses of order Mkk ~ ^/^s^- The 5D SU{Nc) gauge theory 
living on the DA branes has a dimensionful 't Hooft coupling X^d- The dimensionless 
parameter A = As/jMkk has an interpretation as the 't Hooft coupling of a 4D theory 
evaluated at the scale Mkk- In addition to massless gluons, the 4D theory has an infinite 
number of KK modes. Then one can consider the two limits for the coupling at the scale 
Mkk- a — >• and A — t- cxd. If A is small at the KK scale, the KK modes become heavy 
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and decouple from the dynamics at low energies (compared to Mkk)- The only massless 
fields are gluons and fundamental quarks, and the low-energy theory is just QCD, which 
confines at some dynamically- generated scale Aqcd ^ Mkk- The brane system can thus 
be viewed as a particular UV completion of QCD at the scale Mkk- 

On the other hand, one can also consider the opposite limit, with A large at the KK 
scale. Now the Kaluza-Klein tower associated with the X4 direction does not decouple, 
and the theory does not flow to pure QCD in the infrared. However, if both Nc and A 
are large, there is a dual description of the gauge theory in terms of gravity. So at the 
cost of losing a sharp connection to the physics of pure QCD, one can profitably use the 
holographic description to study the strongly-coupled theory, which still turns out to share 
some of the most important features of QCD. This is morally similar to working with 
the strong-coupling limit of a lattice gauge theory, which also allows many simplifications, 
at the expense of losing a sharp connection to the continuum field theory of interest. 
The large Nc and large A limits in the gauge theory are necessary to suppress gs and a' 
corrections in the dual string theory, which then reduces to classical gravity on a weakly- 
curved background. With the optimistic assumption that the A ^ 1 and A » 1 limits of 
the theory are smoothly connected, without any phase transitions at some intermediate 't 
Hooft coupling, it may be hoped that calculations using the holographic theory will give 
results that are qualitatively similar to calculations in QCD"^. With this hope in mind, the 
strong-coupling theory is often referred to as 'holographic QCD' in the literature. 

In the large Nc and A limits, the D4 branes are replaced by their near- horizon geometry, 
with appropriate boundary conditions along for the X4 circle [23]. So long as Nf ^ Nc, the 
D8 branes have no effect on the geometry, and can be treated in the probe approximation. 
In this limit the Type IIA supergravity equations of motion are solved by the metric Gmn, 
dilaton <I>, and RR 3-form C3 fields given in the string frame by [23] 



dU^ 



+ u^dnl 



3/2 /R\3/2 

Uiu) 

^'' = 9s[iy\ m) = l-^, F, = dC,= ^-^e, (2.2) 

where = irgsNclg, rj^y is the flat Minkowski metric, and r = X4 and has circumference 
27r/MKK- Finally, Mkk = 3?7j^^/2i?3/2 ^ ^^^^ -g ^j^^ ^-j^^ element on which has a volume 
V4, and €4 is the volume form on S^. The energy scale of the dual field theory is related 
to ?7M^j^. Modulo some caveats that will be discussed below, U takes values in the range 
(C/kk,oo). The background is topologically R^^^ x I x x [10], with R^'^ x x S"^ 
fibered over the interval /, which is parametrized by U above. 

Flavor fields arise from the D8 and D8 branes[25], which are immersed into the above 
background at the positions t{U — t- 00) = ibJT/2 respectively. There are U{2) fiavor gauge 
fields living on the D8 and D8 branes, which using the AdS/CFT dictionary are associated 
to sources for the currents of the U{2)l x ?7(2)/j chiral symmetry of the dual field theory. 
The embedding function t[U) is determined from the equations of motion of the D8 branes. 



*For some interesting recent developments in this direction at finite temperature, see Ref. [24]. 
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and the D8 branes turn out to connect to each other at ?7 = Ukk- This gives a beautiful 
geometric reahzation of spontaneous chiral symmetry breaking, as discussed in detail in 
Refs. [4, 5]. 

2.2 Region of validity of the supergravity approximation 

To understand the region of validity of the supergravity approximation, some useful rela- 
tions between the stringy parameters Qs, Mkk^Ukk, RJs and the field theory parameters 
A, Nc are 

It turns out that the action depends on Ig only through an overall normalization, and this 
can be used to set Ukk = I/A^kk without loss of generality. Then one obtains the very 
useful relation IM^^P, = X'^. The Type IIA supergravity description is valid provided 
that 1) all curvature invariants are small compared to the string scale, so that a' corrections 
can be neglected, and 2) the effective string coupling, which is controlled by the dilaton, 
remains small, so that string loop corrections can neglected^. In terms of the parameters 
of the dual field theory this translates into the constraint 1/A ^ UMkk *S Nc^^ /X [26]^. 
In this case since the range of U is bounded from below by Ukk, the lower limit above is 
not dangerous so long as A is large, because then a' corrections to the background will be 
negligible. For large enough U, however, type IIA supergravity will not be reliable because 
the effective string coupling becomes large. The fact that the region of validity of Type 
IIA supergravity is energy-dependent can be traced to the field-theory fact that 5D YM 
theory is non-renormalizable. New degrees of freedom, which are not included in Type IIA 
supergravity, become important in the ultraviolet limit. 

In addition to the above issues with the region of validity of type IIA supergravity, in 
order for the holographic description of the theory to be reliable, one must make sure that 
the effective action used to describe the physics of bulk fields living on the fiavor branes 
remains reliable. This latter matter turns out to be a subtle business for baryons, because 
the effective action describing the flavor fields breaks down if the fields vary rapidly on 
distance scales comparable to 1/ a'^/^ . Studies using the fiat-space instanton approximation 
suggest that the size of the cores of baryons in the Sakai-Sugimoto model shrink to the 
string scale, and hence a' corrections to the effective action of the fiavor gauge fields may 
not be under firm theoretical control in the small-C/MKK region; see the discussions in 
Refs. [8-11, 13, 27]. Here we will be concerned with the low-energy properties of baryons, 
which the usual effective field theory ideology implies will only depend on such small- C/ 
issues through a few 'low-energy constants' such as e.g. the mass and moment of inertia 
of the holographic soliton configuration. For the observables we will focus on here, it is 



^Once the dilaton becomes large, one can regain control of the theory by doing a lift to M-theory, so 
that the physics becomes describable in terms of 11-dimensional supergravity 

®Fits to meson-sector data imply that for A^c = 3, A ~ 17, so that a' and gs corrections might be 
important. From a formal perspective it is not obvious why the model works as well as it does in fitting 
the data. 
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expected that all such parameters would cancel out, and we will see that this is indeed the 
case with our approach. 

The treatment of the large U region is much trickier. From the discussion above, 
it is clear that if the holographic description of the theory is to be reliable, one must 
make sure that the observables one is interested in receive no contributions from the region 
C/Mkk > Nt'^/\, where the supergravity theory is not under control. To capture the long- 
range physics of baryons correctly we will see that it turns out to be critical to systematically 
sequester the large U region by introducing a cutoff on U at the intermediate stages of our 
calculations. This brings us to the issue of holographic renormalization. 

In holography, the UV divergences of field theory are traded for IR divergences of 
the on-shell action of the gravity dual due to the infinite extent of the extra holographic 
direction. At the risk of some confusion we will stick with the conventions in the literature 
and refer to such IR holographic divergences as UV divergences. The general setup of 
holography relates the large-C/ behavior of bulk fields to the behavior of sources for various 
operators in the dual field theory. In the well-understood examples of the AdS/CFT 
correspondence, it has been shown that one must impose a cutoff on the holographic 
coordinates at the intermediate stages of calculations in order to systematically remove 
volume divergences in the on-shell action of the holographic theory. Once the theory has 
been regularized and any divergences subtracted by the addition of appropriate boundary 
counter terms, the cutoff can be - and should be - removed. This procedure is known as 
holographic renormalization [28]. 

On general grounds, one expects that the same procedure is necessary for the Witten 
background^. With this motivation, in our calculations we introduce a cutoff on U , so that 
U G {Ukk, Uuv]- However, there is an important subtlety here which is not present in the 
tamer cases to which holographic renormalization is normally applied. In view of the above 
observations on the domain of validity of the supergravity description of D4 branes, C/uv 
must be chosen such that UuvM^k Nc / A so that the behavior of the bulk theory near 
f^max can be trusted. That is, even at the end of a calculation, in principle one should not 
strictly send U^v — ^ co- Instead U^v must be taken larger than the other physical scales in 
the problem, but still remain small compared to Nc^^ /X, which is in principle a large but 
finite number in any given application. Correlation functions in the field theory should then 
be calculated from the on-shell action of the regulated gravity theory, with the boundary 
values of bulk fields at U = U^^ acting as sources for field theory operators. Well-behaved 
observables should have vanishing U^^ dependence as U^vMkk is taken to be large. Indeed, 
it will turn out that working with a fixed cutoff C/uv and only removing the cutoff at the 
last step of the calculations plays a crucial role in the success of the Sakai-Sugimoto model 
in capturing the correct large distance physics of baryons. 

systematic approach to holographic renormalization for non-conformal brane systems (such as the 
Witten background) has been developed only somewhat recently [29-32]. To our knowledge these techniques 
have not yet been applied directly to the Sakai-Sugimoto model. 



-7- 



2.3 Five-dimensional action 

As long as the U{Nf) world-volume gauge fields Am, M = 0, . . . , 8, on the D8 branes vary 
slowly compared to the string scale, their dynamics should be reliably described by a Yang- 
Mills-Chern-Simons action. In the Sakai-Sugimoto model the field theory has an SO{5) 
symmetry, under which the massless fields, which are the same ones that occur in QCD, are 
singlets. We are interested in baryons with the quantum numbers of nucleons, which are 
50(5) singlets, and to this end we follow the approach in the literature of setting A5... ,8 = 
and 95,...,8-^M = 0. It is also convenient to switch to a new holographic coordinate z defined 

by 

= U^k + Ukkz\ (2.4) 

3 /2 1/2 

which takes values in [— ^uvj +-2uv]) with ~ Uuv /^kk' ^^^^ these choices, the Sakai- 
Sugimoto model reduces to a five-dimensional U (Nf) Yang-Mills-Chern-Simons theory in 
a curved background. The only dimensionful scale is MkK) and we choose units such that 
Mkk = 1- The Mkk dependence can always be restored by dimensional analysis. The 
action is then[4, 5] 

/2uv 
d'^xdz Tr 

Here k = AA''c/(2167r^), and the functions h{z) = (1 + z"^)^^^'^ and k{z) = 1 + encode 
the metric of the curved background. The field strength is given by = dA — lA A A 
and UJ5 = Tr {AT'^ - iA^T/2 - A^/IO) is the Chern-S imons five-form. In particular, we 
consider the case that Nf = 2. 

Finally, Sb stands for a set of boundary terms at 2; = Zuv, which must be constructed 
from the boundary data in such a way that the total action is gauge and diffeomorphism 
invariant [28] . The coefficients of some of these terms may be fixed by the requirement that 
the on-shell action evaluated at while the coefficients of 

others may not be fixed a priori. The latter sort of terms are called 'finite counterterms' in 
the literature, and the freedom to choose their coefficients reflects the freedom to choose 
a renormalization scheme. Physical observables must of course be scheme-independent, 
but the detailed connection between 'bare' model parameters and physical observables 
is in general scheme-dependent. The observable we will be calculating, Eq. (1.1), is a 
pure number, and by construction cannot depend on any model parameters in any theory 
consistently implementing the expectations of chiral perturbation theory. Thus it also 
cannot depend on any of the boundary counterterms in S^, and as a result in this paper 
we will not need to work with Sb in an explicit form. It is nonetheless important to keep 
in mind that in principle Sb / 0, and its detailed form should be expected to affect the 
mapping between the model parameters and field theory observables. 

Ultimately we seek to compute the electromagnetic form factors of the baryon, which 
are evaluated as matrix elements of vector and axial isospin currents. The currents can be 
computed using the standard ideas of gauge-gravity duality. The gauge theory lives on the 
UV boundary (z = z^v in our case) of the gravity theory. Each gauge-invariant operator 



1 



/i(2) + kiz)J=-; 



+ 



247r2 



W5 + Sb- 



(2.5) 
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in the gauge theory is associated with a dynamical field in the gravity theory, with the 
boundary value of the bulk field acting as a source for the operator. The partition function 
of the field theory is equated to the partition function of the gravity theory with all bulk 
fields evaluated on-shell. In particular, since the source for a conserved current in a field 
theory is a gauge field, the bulk fields associated with conserved currents are gauge fields, 
which are precisely the A fields in the action above. As shown in e.g. Ref. [13] in the 
Sakai-Sugimoto model the isovector Jy, axial isovector J a, isoscalar Js and axial isoscalar 
Ja currents are 

(2.6) 
(2.7) 
(2.8) 
(2.9) 

where the bulk fields are evaluated on solutions to the bulk equations of motion, and we 
split A (and hence J^) into SU{2) and U{1) pieces as A = A + \'^2A. 

2.4 Baryons as solitons 

Baryons are encoded in soliton solutions of the model. In the D-brane construction, a D4- 
brane wrapped on S"^ provides a baryon vertex in the bulk [7], and appears as instanton-like 
solitonic configurations in the gauge theory on flavor D8-branes[4, 5, 8-10, 33]. In the effec- 
tive 5D YM+CS theory, the baryon solitons carry unit instanton number on xi,X2, X3, z ^ . 
We start from the static solutions for simplicity. Collective coordinate quantization of the 
soliton, which is necessary to pick out the physics of protons from the physics of solitons, 
will be considered later. 

There are a number of symmetries we can take advantage of to aid in finding the 
solutions. On general grounds, we expect that the ground-state baryons in the exact 
isospin limit will be spherically symmetric. On the gravity side, the action has an 5*0(3) 
spatial symmetry, and so it is natural to expect that the minimum-energy static soliton 
will be SO{2>) symmetric, provided that the boundary conditions also have this symmetry. 
In the usual approach to baryons in the Sakai-Sugimoto model, which relies on the flat- 
space instanton approximation, one essentially seeks an SO {A) symmetric configuration, 
since there are some arguments that an approximate 50(4) symmetry appears at large 
A. Here we want to avoid relying on the flat-space instanton approach, and do not ask 
for 50(4) symmetry. In any case, the 50 (3)-symmetric ansatz is the only self-consistent 
ansatz one can impose a priori (without making approximations), since once we add the 
point at spatial infinity the topology of our manifold is 5^ x /, not 5^. On top of this any 
putative 50(4) symmetry of the geometry is broken by the warp functions h{z), k{z). 



* For some recent discussions of 5D instantons in supersymmetric YM+CS theories in flat space, see 
Refs. [34, 35]. 
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The correct procedure for imposing spacetime symmetries on gauge fields was elu- 
cidated in the work of Forgacs and Manton[36], following Witten[37]. In the context of 
holographic baryons, this procedure was first applied by Panico, Pomarol and Wulzer[20- 
22], whose notation we will follow. In particular, it turns out that consistently imposing 
an S0{3)s rotation symmetry on a U{2) gauge field Am necessitates demanding that Am 
be invariant under the diagonal combination S0{3)s + 50(3)^ transformations, where the 
second factor is the 5*0(3) subgroup of the gauge group U(2). This procedure produces a 
self-consistent S0{3) symmetric ansatz for Am that can be plugged back into the action. 

For time-independent configurations, the YM-CS action has a Z2 symmetry which 
implies that we can set Aq and the spatial components of A to zero[9, 21]. The 5*0(3)- 
symmetric time-independent ansatz is then given by[20-22, 36, 37] 

At = A, — , Ao = s. (2.10) 
r 

where (j)2, A^, Ar, s are all functions of r and z; i, j = 1, 2, 3 run over the spatial indices, 
and a = 1,2,3 is an isospin index. The fields </>i,(?5>2 can be naturally packaged into a 
complex scalar field (j) = (pi + i(p2 which has unit charge under the Abelian gauge field 
with components Az,Ar, so that -0^0 = + iAf^(p, where now fj, = z,r. The reason for 
the presence of the 2D Abelian gauge field in the ansatz can be traced to an unbroken 
Abelian subgroup of original SU{2) gauge symmetry in the symmetry-reduced ansatz. 
Plugging in the ansatz, the mass-energy of the reduced system can be written as 

M = Mym + Mcs + Mb, (2.11) 

where / dtM = 850,! dtMf, = Sb, 



00 



Mym = 167rK dr dz 

'0 J~Zuv 



h{z)\Drcl)\^ + k{z)\D,^\^ + \r^k{z)Fl, (2.12) 



+ ^^h{z){l - \cp\''f - y {h{z){drsf + k{z){d.sf) 



(2.13) 



and 



/>00 /-Zuv 

Mcs = 167rK7/ dr dz s ef"" [df,{-icl)* D^c/) + h.c) + F,,^], (2.14) 

with 7 = Nc/{IQ'7t'^k) = 277r/(2A). The expression for the topological charge in terms of 
the reduced ansatz is 

Q=^J drdz (e^'^'d^ [-icP*D^cP + h.c] + e^^F^,^). (2.15) 

Since by definition Sb is a pure boundary term, it does not contribute to the bulk equa- 
tions of motion. We will also shortly see that with our choice of boundary conditions Sb 
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manifestly cannot contribute to Q. The equations of motion that fohow from extremizing 
Mym + Mcs are 

= Dr {h{z)Dr(t>) + {k{z)D,<p) + ^(^(1 - 101') - iie^^d^sD^cf), (2.16) 
{) = dr {r^k{z)Fr,) - k{z) {i(l)*D,(j) + h.c.) - -fe'''drs{l - |(/>p), (2.17) 
= d, {r^k{z)F,r) - Hz) ii^*Dr<j) + h.c.) - 7e^"a,s(l " I0P), (2-18) 
= dr {hiz)r^drs) + d, {kiz)r'^d,s) + 7e^^ [d^{-i^*D^^ + h.c) + F^,^] . (2.19) 

The equations of motion are a set of coupled second-order non-linear PDEs, and a 
general analytic solution of them seems out of reach. In general one would have to construct 
the solutions numerically. Fortunately, given our goal of investigating the long-distance 
properties of the baryon, we only need to know the behavior of the solution at large 
distance r. For this we only need to solve the equations of motion to the first few orders 
in a 1/r expansion, and this can be done analytically. Implicit in our analysis will be the 
assumption that there exists a global unique and well-behaved solution to the boundary 
value problem with Q = I. 

We wish to choose boundary conditions which will give finite-energy solutions with 
charge Q = 1. The integrand in the Eq. (2.15) is a total derivative, and Q can be written 
as 



Q 



1 

2^ 
1 

2^ 



dr 



dz 



-i(t)Dr4>* + i(p*Dr<f>) + Ar 



{-i4>D-A* + i(p*D,cP) + A, 



-H~-2uv 
- ^uv 

r=R 
r=0 ' 



(2.20) 
(2.21) 



with R 



oo. 



To find the large-r solutions, we find it convenient to work in the 2D 



Lorentz gauge 9^^!^ = 0, and choose the boundary conditions for the fields 4),A^,s in 
such a way that the only non-trivial contribution to Q comes from the boundary at r = 0. 
These boundary conditions are summarized in Table 3.2. With this choice of boundary 
conditions Sb certainly cannot contribute to Q. The 5D SU{2) gauge fields are singular at 
r = in this gauge, but of course gauge-invariants like the field strength remain smooth 
at the origin. 



r — )• oo 


r = 


Z — ibZuv 


01 = 


01= sin (2^) 


01 =0 


02 = -1 


02 = -COs(^) 


02 = -1 


A, = 


A — ^ 


dzA, = 


drAr = 


drAr = 


Ar = 


S = 


S = 


s = 



Table 1. Boundary conditions for (p — (pi +i02, Afj_ and s in the 2D Lorentz gauge = 0. Here 

L = 2zuv is the extent of the holographic direction with the UV cutoff Zuv 
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3 Large r properties of baryons 



We are interested in the solutions at large r. If r is much larger than the characteristic 
physical scale of the model, M^^, the equations of motion should linearize. This motivates 
looking for asymptotic solutions to the equations of motion as a Taylor series in 1/r, with 
coefficient functions that depend on z. If the coefficient functions of the Taylor expansion 
can be solved for order by order, the linearization of the equations of motion at large r will 
have been demonstrated self-consistently. It is important to note the physics behind these 
statements. Seeking solutions to the equations of motion for a soliton which are Taylor 
series in 1 /r at large r is tantamount to asking that the profile function of the soliton have 
a long-range pion tail. Showing that the coefficient functions of the Taylor expansion are 
non-trivial and can be solved for order by order then shows self-consistently that the profile 
function indeed has a long-range pion tail, and hence implies that the baryon is surrounded 
by a pion cloud. 

3.1 Static solutions 

To set up the expansion we write (pi, (1)2, A^, Ar, s as 

oo ^ oo ^ 

0i(r,z) = ^4")(z)-, </>2(r,.) = -l + 5^4")(z)-, (3.1) 

n=l n=l 
oo ^ oo ^ oo ^ 

AAr,z) = Y,A(rHz)-, Ar{r,z) = Y,Ai-Hz)-, s{r, z) = Y,s^^\z)-. (3.2) 

n=l n=l n=l 

This ansatz is chosen to automatically satisfy the boundary conditions at r = oo. The 
equations of motion then reduce to a set of coupled ODEs for the coefficient functions 
(f)^^\z), 02"''(2;), ^i"^(-z), A^^\z), s(")(z). It is then possible to solve for the coefficient func- 
tions order by order in 1/r using the boundary conditions at z = izuv Counting boundary 
conditions we see that we have enough data to fix all of the solutions order by order in 
terms of a single undetermined constant of integration, which we call /3 below. 

The role of /3 is to parametrize the dependence of the large-r solution on the rest of 
the solution. To see this, note that the physical boundary value problem that determines 
the full soliton field configuration depends on two parameters, z-^v and 7 ^. If one makes 
the standard physical assumption that there exists a unique well-behaved solution to the 
boundary value problem with a given and 7; then will be fixed by matching to full 
solution, and will be a function of Zuv and 7. 

Solving the equations of motion at large-r order-by-order in power series in 1/r, we 
find that the first few terms in the expansion are 

, ^ ^ tan -^(zuv) ' /„ o\ 

(pi = 5 [6.6) 



(3 I -3 (-1 + z"^) Zuv tan"^ z + z [{-6 + z"^ + 2zl^) tan"^ Zuv + 32:uv log 



l+z'^ 



tan -'-[zuv]?''' 



^Assuming that the topological charge is fixed to Q = 1. 
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-1 + 



tan 1 [zuv] 



(3.4) 



p P (dzz^y tan-1 [z] + (3 - 3z^ - 2zl,) tan"! [z^y] - 3z^y (l + log [^] ) ) 



-2zP + 2^uv/3tan-i[,] 
A _ ^ tan ^[zuvj , 

-fir — n r 



tan -^[^;uv]?'^ 



(3.5) 
(3.6) 



4P (-3 (-1 + z^) Zuvtan-i[2] + z ((-3 + z^ + 2zl,) t&n-^[zuy] + 3zuvlog [rf^])) 



tan~^[zuv]r^ 

P^7zly{tan-^{z)^ - 6tan-i(z)2tan-i(^uv)^ + 5tan-^(z)tan-^(zu 

2tan-i(zuv)^r9 



(3.7) 



It is crucial to work with a finite cutoff z^y to obtain these solutions. If one tries to take 
Zuv ~^ oo when looking for the solutions, only the trivial solution to the boundary-value 
problem will be found. To see this from the solutions above, consider for instance (pi. At 
large r, 

In the large z^y limit, $ — >■ — ^^^tan(2;), which solves the differential equation for (f)i. 
However, the only way to satisfy the boundary condition c/ii — )• as z — > oo is to set /3 = 0. 
So if one tries to work without a cutoff on z, as has thus far been the approach in the 
literature, the only solutions for the boundary value problem which are power series in 1/r 
that can be obtained are the trivial ones. One then cannot capture all of the long-distance 
physics. With a cutoff, however, one obtains non-trivial solutions to the boundary value 
problem order by order in a 1/r expansion, and the model has a chance of capturing the 
physics of interest. 

With the above solutions in hand, we can evaluate the asymptotic on-shell action 
density for static solutions at large r, with the result that 

(. - ..)i_ ^ / ,uz { + on/r^)] = / { + ayr^)] 

(3.9) 

^ /^*{ ^^^1^ + ^(1 A')}' (3-10) 

where in the last line we took z^y S> 1. There are then two possibilities for the behavior of 
the on-shell action as z^y oo: either P behaves as /3 ~ l/^^uv at large z^y in the full Q = 1 
finite-energy solution, or it does not. If /3 behaves as /3 ~ l/^uv^ then the on-shell action 
will have no large z^y divergences at large r, and Sb can only contain finite counterterms. 
Otherwise, the on-shell action without ^5 may be divergent at large z^y, and then Sh would 
also contain terms that diverge at large z-ay in such a way that the full on-shell action is 
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finite. A very interesting task we leave to future work is to work out the beliavior of the 
full solution and understand which of these options is the relevant one by solving the full 
PDE and examining the asymptotic behavior of the solution. Fortunately for the current 
analysis, the evaluation of the form- factor ratio is independent of these interesting and 
subtle issues. 

3.2 Rotating solutions 

To extract information about states with the quantum numbers of the proton, we must 
follow the procedure of collective-coordinate quantization [38]. This is because in general, 
single-particle states associated with solitons are described as quantized time-dependent 
fluctuations around the static soliton solution in the zero mode directions. To perform 
collective coordinate quantization, one allows the soliton to slowly rotate in the zero-mode 
directions. Then the gauge field components Aq and Ai,Az must be turned on so as to 
parametrize the collective motion of the rotating soliton. In the limit of very slow rotation 
(ultimately justified by the large Nc limit, which makes the moment of inertia of the soliton 
large), the shape of the soliton cannot be affected by the rotation, so we expect the soliton 
to keep its SO{3) rotational symmetry. This S0{3) symmetry can then be used to constrain 
the form of Aq and Ai,Az. The general SO{3) symmetric ansatz appropriate to a soliton 
rotating with a constant angular velocity k is given by [22] 

+ v{k-x)x'', (3.11) 

Ai = ^(ki-ik- x)xi) + BJk ■ x)xi + Qeibck^x"", (3.12) 
r 

A, = B,{k-x). (3.13) 

The way the new fields xi,2, P, Br^z,Q,v come into the symmetry-reduced action is con- 
strained by the residual gauge symmetry. The field X = Xi + transforms the same way 
as 4>, and couples to the gauge field A^. The time-dependent ansatz has a new residual 
C/(l) gauge symmetry, associated with chiral C/(1)l,h transformations of the form = g, 
QL = with 

g = exp ^i^{r, z){k ■ x)J . (3-14) 

Under this transformation only and p transform non-trivially, with transformation rules 

B^^B^ + d^C, p^p + C. (3.15) 

Evaluated on the 50(3) symmetric time-dependent ansatz, the YM-CS Lagrangian 
becomes 

C = -M +^kak'' (3.16) 



^0 



+ X2 



[x^x^ 
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to first order in k, where M is the mass of the static configuration and was given in 
Eq. (2.11), and the moment of inertia A can be written as 



A 



dr 



dz 



h{z){DrpY - k{z){D,pY - r'h{z){drQY - r'k{z){d,QY 



2/i(z)Q2 



h{z)Bl, + r''hiz)\DrX\^ + r^kiz)\D,x\^ + -h{z)idrvY + -kizYd^vf 



+ h{z){\x\^ + v^){l + \<P\^) - Av^pXp + 7( - 2e^'''D^pXp{D,cl))p 

+ 2e^^5^(rQ)xpeP^(Z),0), - v{e>'''B^,m^ - l)/2 + rQe'^^ ^ + 2rQe^'' D^pd^s 



(3.17) 

where p,q = 1, 2. Note that A ~ Nc- The form of Eq. (3.16) is that of a Lagrangian for a 
rigid rotor of mass M and moment of inertia A, with collective coordinates ka, a = 0, 1, 2, 3. 
The equations of motion that follow from this Lagrangian are 



h{z) 



h{z) 



l)Brz + 2rQFr,] = 0, 



^Drir^Drx) + j;^MHz)D,x) + \{2vcl> - (1 + |0|2)x) - 

\dr{r^drQ) + -^a,(A;(z)5,g) - 4q 



e'^"(L»^0(ia,Q + L>,p) = 0, 



2r/i(z 



/i(z) 

-e^'^'iiiD^ct^iD.x)^ + h.c.) + F^,(2t; - x4>^ - 4>X^)l2 - 2D^pd,s] = 0, 



7 



9.(A.P) + ^9.(.(.)Z).P) ^^^^^ 



'^'^[(I)^,^(Z),X)^ + /i-c.) + iF^u{<Px^ - X0^)/2 + 29^(rQ)5,s] = 0, 



^ a,(fc(z)S,,) + \DrP + ^T^[((X - ^^'A)(^.'A)^ + h.c.) + (1 - |</>|2)5,t; - 2rQd,s\ = 0, 



/i(z 



r'^h{z) 

^^[((X - ^;(/')(I)r</')^ + h.c.) + (1 - |(/<|2)a,7; - 2rQ9,s] = 0. 



(3.18) 



The boundary conditions for the new fields are given in Table 3.2. 



r 


r = oo 




Xi = -sin(^) 


XI = 


XI = 


X2 = cos(^) 


X2 = 1 


X2 = -1 


= 


5, = 


= 


drBr = 


drBr = 


Br = 


/5 = 


p = 


,9 = 


u = -1 


t; = -1 


v = -l 


Q = 


Q = 


Q = 



Table 2. Boundary conditions for X — Xi+ iX2; -B^, Q, and p in the 2D Lorentz gauge d^B^^ — 0. 



Using the same approach as in the static case, it is straightforward to solve the equa- 
tions of motion (3.18) to the first few orders in a 1/r expansion. The solution depends 
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on the matching parameter (5 already seen in the static solutions, and on an additional 
matching parameter /3. The results are 



_ . o z,.,/9tan-i(^) 
Xl = 5 ^ 



(3.19) 



P (-3 (-1 + 2^^tan-i(z) + z (-3 + + 2^2^) tan-i(zuv) + 32;Zuvlog [r^]) 



tan ^(zuv)''^ 



X2 = l + 



^ tan~-'(2uv) 
-3(-l+z^)z„vtan-^{z) 



+ 



1+^2 



g ( n::-"::^' + (3 - S''^ - ^^v) - ij^, (i + log 



= -l + 0(l/r^) 

_ 4v/3^7 (tan-i (z)^ - 6 tan"! {zf tan"! (zuv)^ + 5 tan'^ (zuv)^) 
^~ 2tan-i(zuv)V^ 



1+2^ 



1 + ^,1 



tan H-^uv) 



(3.20) 
(3.21) 



(3.22) 
(3.23) 

(3.24) 



P (^-3 (-1 + 2^) tan-i(z) + z (-3 + + 22;2^) tan-i(zuv) + 3zZuv log 

tan~-^(2uv)^^ 



Evaluating the asymptotic on-shell moment of inertia at large r, we obtain 

3 7 I «;(2;)(tan -^^uv)^ A;(2;)(tan Zuv) I 

^ IGttk f 2(2uv/3)2 6(2:uv^)2 6^ 

3 I tan-i(^uv)r2 ^ tan-i(^uv)r4 J 

^ ^ + '-"^ + 0(1A^)| . (3.25) 

3 I Trr"' 7rr* I 

We observe that for the rotating solution to have finite energy, either /3 must go to zero at 
least as fast as ~ 1/zuv, or the /3 contribution must be cancelled by a contribution from 
the boundary action Si,. Fortunately, as will be seen shortly $ does not contribute to the 
leading large-r behavior of the form factors, and so we can defer an investigation of the 
behavior of $ to future work. 
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3.3 Evaluation of the form factors 

In Ref. [15] it was argued that X"PT constrains the Fourier transforms of the electromag- 
netic form factors of the proton to take the large-r forms 

Here 7,^ ~ is the pion decay constant, gA ~ is the axial couphng constant, and 
A = (Ma — Mn)/ Aqcd ~ l/-^c is the nucleon-delta mass sphtting. The ratio in Eq. (1.1) 
is constructed to make the low-energy constants cancel, so that it should take the same nu- 
merical value in any theory with baryons and the same pattern of chiral symmetry breaking 
as in QCD. Of course, this assumes that a particular model in question does not suffer from 
the issues discussed in the introduction to do with the an accidental (from the perspec- 
tive of x-PT) suppression of some of these low-energy constants. If such suppressions are 
present and come in an unfortunate pattern, the long-range physics might end up being 
qualitatively different. 

The quantization of the soliton system is standard[38], and is discussed at length in 
this notation in Ref. [22]. For our purposes, the important outputs of this analysis are the 
following identities for the matrix elements of the collective coordinates on the subspace of 
nucleon states: 

( T^- Ua'^U^a") = -^S'^P, (3.30) 

( Tr Ua^Xb{k ■ xt/V") = (3.31) 

(3.32) 

where 3°" and 1°" are the expectation values of spin and isospin respectively. 

Armed with these identities, we can plug the expressions for the isoscalar and isovector 
currents in Eq. (2.6) into the definitions of the position-space electromagnetic form factors 
in Eq. (1.2), and thence obtain simple expressions for the electromagnetic form factors in 
terms of the fields parametrizing the symmetry-reduced ansatz in Eqs. (2.10), (3.11). The 
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result of these manipulations is 



&E\r) = -i-^ ' (3-33) 

Gil\r) = [rHz)d,Q 1^ > (3-34) 

G'E=\r) = [k{z){d,v - 2{d,X2 - ^.Xi))]X ' (3-35) 
4 

--/ 
9 



Gil\r) = --K [k{z){d,ct>2 - A,h)]X^ , (3.36) 



Now we are basically done, since we have everything we need to check the model- 
independent relation Eq. (1.1). Plugging in the large r solutions in Eqs. (3.3), (3.19) into 
Eq. (3.33), and extracting the leading terms at large r and z^v, we get 



Gi7°(r) 
&E=Hr) 
Gitir) 



^7=0 ^ 32Atzg,/33^ 



_16kz^ 
97rr4 



Note that in all of these expressions /3 enters in the combination f3z^^. The form factors 
are observables, and so cannot depend on a 'UV cutoff or any choice of renormalization 
scheme in the gravity dual. However, the detailed expressions of physical observables in 
terms of model parameters can certainly depend on the choice of scheme. There are two 
generic ways these expressions may become independent from the cutoff. The first is if /3, 
which is determined by the global solution to the boundary value problem for the Q = 1 
soliton, scales as /3 ~ l/^uv when z^^ ^ I. Alternatively, f3 may scale differently, and there 
may be contributions from Sh to these relations that we did not explicitly write. To work 
this out completely, one could for instance verify that when Eqs. (3.37) are rewritten in 
terms of /jr, ffA, and A, all scheme and cutoff dependence disappears, so long as fn,gA, and 
A are themselves calculated using the same cutoff we used above. 

Fortunately, for our main goal of calculating Eq. (1.1), these difficult issues do not 
have to be resolved. Assembling the ratio, we find that 

hm r^^Jf-^ = 18, (3.38) 

in the Sakai-Sugimoto model, as advertised. As we have been foreshadowing, all of the 
model parameters and the associated subtleties cancel from the ratio Eq. (3.38), as they 
must in any model that consistently implements the structure of baryon-pion physics ex- 
pected from large Nc chiral perturbation theory. The reassuring implication is that the 
structure of the long-range physics of baryons in the Sakai-Sugimoto model is the same as 
in large Nc QCD. 
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4 Discussion 



Our results provide strong evidence that baryons in the Sakai-Sugimoto model obey the 
model-independent relation Eq. (1.1). To see this, we focused on the behavior of the 
soliton solution of the model at large distances r. Armed with asymptotic solutions of the 
equations of motion obtained in power series expansion in 1/r at large r, we evaluated the 
large r contributions to the electromagnetic form factors of the proton. The form factors 
behaved in precisely the way expected from X"PT, and as a result the large r form factors 
obeyed Eq. (1.1). So the issues of the Sakai-Sugimoto model appearing to fail to meet the 
expectations of X"PT raised in Ref. [15] seem to be connected to subtle issues with the 
applicability of the standard approach in the literature of using the flat-space instantons 
approximation, rather than a deep problem with the model itself. 

It remains a very interesting task to understand exactly why the standard flat-space 
instanton approach gives misleading results for the observables discussed here, and to 
clarify the circumstances in which the standard approach should be reliable, since it leads 
to apparently reasonable predictions in many cases. Here we note a couple of tentative 
possible reasons for the discrepancies between our results and those in the literature. The 
first one has to do with the standard approach in computations in the Sakai-Sugimoto 
model, where z^^ is taken to infinity from the start of the calculation. As we saw in the 
direct construction of the large r solutions, sending Zuv — ^ oo too early can cause one to 
misidentify the asymptotic structure of the soliton fields. Another possible issue concerns 
the usual approach for extracting the behavior of the currents from the flat-space instanton 
solution. For this it is necessary to know the behavior of the fields at large z, but there 
the fiat-space instanton approximation cannot be trusted. The standard approach is to 
argue that the soliton solution linearizes already at z <^ 1, and the flat-space instanton 
solution is matched to the solutions of the linearized equation of motion for z ^ 1. But as 
we saw in the direct solution at large r, while the system certainly linearizes in the sense 
that an order-by-order solution is possible, the non-linear terms in the full equations of 
motion play an important role in determining the solutions. This suggests that the task of 
connecting the flat-space instanton to the behavior of the soliton fields at large z may be 
rather subtle. 

There are many important directions for future work. Perhaps the most urgent task 
is to find a numerical Q = 1 solution to the equations of motion, and match it both to 
our large r solution and to the r,z A~^/^ flat-space instanton solutions in the literature. 
Among other things, this would greatly help in understanding the connections between the 
approach taken here and the flat-space instanton approximation. 

Another urgent and possibly related task is to undertake a systematic study of holo- 
graphic renormalization for the Sakai-Sugimoto model. Thanks to the special nature of the 
observable we were aiming at [Eq. (1.1)], we were able to get away with being somewhat 
cavalier in our treatment of such issues in our current work, but the problem clearly cries 
out for a careful treatment. A detailed understanding of holographic renormalization and 
of the boundary counterterm action Sb may be important both for constructing numerical 
soliton solutions, and for understanding the domain of validity and proper interpretation 
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of the previous calculations in the literature, both in the meson and baryon sectors. 

Finally, there may be some phenomenological implications of our results. It would 
be interesting to revisit the many baryonic observables already computed in the Sakai- 
Sugimoto model, and see whether going beyond the flat-space instanton approximation 
improves the phenomenological match to real-world data. It is possible that such an ap- 
proach may shed some light on the large-distance interactions between baryons in the 
Sakai-Sugimoto model, which are known to be difficult to accurately capture with the flat- 
space instanton approximation[27, 39-44]. As another direction, it may be interesting to 
consider the implications of our results for the description of baryons with valence strange 
quarks away from the SU (3) flavor limit, which are most generally approached using the 
bound-state approach [45, 46] in the context of the Skyrme model. A trial study in this 
direction in the Sakai-Sugimoto model has been performed in Ref. [47], where flat-space 
instanton approximation was adopted. Since it seems that contributions of pion clouds 
may be significant for the bound-state approach, it will be interesting to reexamine these 
calculations in light of the present results. 
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